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ABSTRAC T 


The purpose of this thesis was to postulate a method of long-range 
inertial navigation, using only acceleration inputs, and to se Bening theo- 
retically the response of this system to external disturbances, using vari- 
ous types of mechanizations. The problem was divided into one of track 
control and of range indication, using a great-circular path between the 
points of departure and destination. 

The response of the track control system was examined with mecha- | 
nization equations of as high an order as the fifth, and was found to improve 
as the order of the mechanization equation increased. Response to impulse 
type wind acceleration disturbances was satisfactory, but long-period sinu- 
soidally varying winds caused excessive errors. It appears probable that 
additional feedback loops, in conjunction with a mechanization equation of 
the fifth or higher order, will solve the track control problem satisfac- 
torily. 

In the range indication system, it was found to be impossible to re- 
move any forcing function terms of higher order than the time rate of 
change of acceleration. It appears probable that the range indication sys- 
tem, as postulated, will prove to be satisfactory, especially if the sensi- 


tivities are made variable, as a function of the angle of the input pendulum. 
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CHAPTER I 


INTRODUCTION 


The importance of developing a long-range self-contained system of 
automatic, mid-course, navigation for the control of guided missiles can- 
not be overestimated (Fig I-1). As long range flights of piloted aircraft 
over enemy territory become increasingly hazardous, because of the neces- 
sity of travelling without fighter protection toward heavily defended cities 
provided with elaborate warning nets, automatic navigation can provide a 
method of delivering bombs without risking the lives of highly trained 
pilots. 

With the development of long-range jet-propelled missiles, automatic 
navigation should eventually permit the United States to bomb any place in 
the world from bases located within its territorial limits. 

For such a system of automatic navigation to be successful, it must 
have a degree of accuracy comparable to that achieved by human pilots. It 
must not be easily susceptible to "Jamming" by enemy action. It must op- 
erate successfully at supersonic velocities, and in the upper reaches of the 
earth's atmosphere. 

At present, the problems of automatic long-range navigation are under 
consideration by several groups. Among these are: 


Baird Associates, Inc. 

Hughes Aircraft Co. 
Instrumentation Laboratory, M.I.T. 
Kolisman Instrument Company 
North American Aviation, Inc. 
Northrup Aircraft, Inc. 
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7. Rand Corporation 


8. U.S. Army Research and Development 
Sub-office (rocket), Fort Bliss, Texas 


Little work, however, appears to have been done on the derivation of 
theoretically optimum control equations or operating parameters. For 
this reason the writers determined to examine a long-range automatic 
navigation system now under development at the Instrumentation Labora- 
tory, M.I.T., as a means for exploring theoretically some of the general 
problems of automatic navigation. The Instrumentation Laboratory proj- 
ect, which is being executed under USAF Contract W33-038ac-13969, is 
designated "An Automatic Navigation System - Project Febe," and will be 
the subject of a report soon to be published. This project has for its pri- 


mary purposes to determine: 


1. the feasibility of long-range automatic navigational 
guidance of bomber airplanes 


2. useful design parameters for a serviceable system 
for military use. 


Considerations other than military and tactical indicated that the sys- 
tem should employ solar tracking, a magnetic azimuth system, and a con- 
stant ground speed. The system is allowed to make no contact with the 
earth. Altitude is determined by the use of a barometric altimeter; the 
only inputs to the system are celestial observation, the magnetic field of 
the earth, and the various accelerations experienced by the airplane. 

This system, which has been installed in an Air Force B-29 bomber, 

' causes the airplane to fly a definitely programmed great circle course at 
constant ground speed. Vertical accelerations are assumed to be of such 
small importance that they can be neglected. Figure I-2 shows a functional 
diagram of the Febe automatic navigational system. During the period from 
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June to August of 1948, the writers of this thesis were concerned with the 
Febe project in an under-instruction status. Through the courtesy of Dr. 
C.S. Draper, they were also permitted to attend a Seminar on Automatic 
(Celestial and Inertial) Long-Range Guidance Systems conducted at the 
Massachusetts Institute of Technology by the Scientific Advisory Board 
from February 1 to February 3, 1949. 

This background naturally led the writers to concern themselves with 
the form which a practicable service guidance system might take. The 
Febe system presents several technical and military problems, which 
arise principally in connection with the following: 

1. Constant ground speed 
2. The necessity of solar tracking 
3. Magnetic azimuth input 
4. Weight and size of components. 

Since a successful automatic navigational system is ultimately des- 
tined to fly in a long-range high-speed missile, it must be assumed that 
the weapon will be operating at close to the maximum range permitted by 
its size and fuel capacity. The navigational system, therefore, should not 
seriously lower the fuel economy. This suggests a system which will fly 
at constant, or nearly constant, airspeed. Furthermore, if the power plant 
consists of an athodyd, the airspeed must be maintained very nearly con- 
stant by the fundamental limitations of this type of propulsion. 

The limitation of the Febe system to a constant ground speed was 
largely dictated by the requirements of the azimuth system. The magnetic 
input could be eliminated through the use of two an hiaaliik ts telescopes. 
The problems of celestial tracking become increasingly severe, however, 


as missile speed is increased. At supersonic and near-sonic speeds, it is 
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no longer feasible to have an astrodome, in which to house the tracking 
unit, projecting from the fuselage. Thus, tracking through a flat window, 
with the attendant difficulties which this entails, becomes a necessity. 
Also, at high speeds, thick boundary layers and intense heating will exist 
along the surface of the missile, greatly complicating optical problems. 
Finally, the inclusion of a celestial tracking system materially increases 
the space and weight required by the navigational system — and weight and 
space factors become increasingly important with increasing range and 
bomb load of an aircraft. It appears, however, that approximately a thirty- 
fold decrease in the uncertainty levels of existing gyros would permit the 
maintenance of an inertial coordinate reference system within the missile 
through the use of such gyros alone, thus eliminating the use of celestial 
tracking. 

Great progress is being made in the improvement of gyros. The work, 
for example, that is being done by the group at the Massachusetts Institute 
of Technology under Dr. C. S. Draper shows promise of obtaining the nec- 
essary accuracy within the next few years. 

The model of the Febe automatic navigational system presently in- 
stalled in a B-29 is very bulky, and weights 1917 lbs. (See Fig 1-3). 

It practically fills the after pressurized compartment of the airplane. Of 
course, it must be remembered that no effort was made to decrease the 
size of this equipment, or to keep its weight at a minimum. Nevertheless, 
a missile held to a constant airspeed, rather than to a programmed ground 
spped, could, by using gimbal solvers, greatly simplify the elaborate trigo- 
nometric and speedcomputers which Febe requires. Much of the work of 


the computation systems might, of course, be performed on the ground 
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prior to the start of the flight, and be fed to the system through tapes, but 
this would introduce further complications. 

From all of these considerations, a self-contained inertial-gravita- 
tional system flying a great circle course at nearly constant airspeed 
appears to provide one of the most satisfactory solutions, if such a system 
can be implemented. 

After this examination of the Febe system, there remain, then, for 
use in the purely inertial system, a gyro-stabilized inertial platform, a 
clock mechanism to remove the earth's diurnal rotation, and two mutually 
perpendicular single-degree-of-freedom pendulous accelerometer units. 
In addition, there are gimbals and servomechanisms isolating the inertial 
platform from the motions of the missile containing it, and orienting the 
platform so that its axis is parallel to the polar axis of the earth, so that 
a reference vertical parallel to that of the point of departure (or any other 
convenient reference direction) is maintained independent of the motion of 
the missile. (See Figs I-4, I-5 and I-6.) 

The three integrating gyros which are used to maintain a gimbal- 
mounted inertial platform "fixed" in inertial space, parallel to the earth's 
equatorial plane, are mounted so that they detect motion of the platform 
with respect to inertial space about the polar axis of the earth, and about 
two mutually perpendicular axes parallel to the plane of the equator (Figs. 
I-7, 1-8). Servomechanisms, controlled by these gyros, rotate the gimbals 
so that the inertial platform remains fixed in inertial space. The earth 
platform, mounted on the inertial platform with the identical polar axis, is 
rotated about that axis once in twenty-four hours, to transform the inertial 


reference system of the inertial platform into an earth reference system. 
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Upon the earth platform, in turn, the controlled member is mounted 
in gimbals (Fig I-9). The axis of the outer gimbal of this member is 
made to lie in the plane of the great circle connecting the points of de- 
parture and destination on the surface of the earth, considered as a 
sphere.* Then, about the inner controlled-member gimbal axis, the con- 
trolled-member is moved until it lies parallel to the surface of the earth 
at the point of departure. The inner controlled-member axis is placed 
perpendicular to the programmed great circle path so that, as the missile 
moves along the programmed path, the controlled member, by rotation 
about the inner axis alone, can be maintained at all times parallel to the 
surface of the earth below the missile. The angle through which the con- 
trolled member has rotated will then be a direct reading of the distance 
travelled by the missile along the programmed great circle. At some pre- 
determined angle between the controlled member and the gimbal connect- 
ing it to the earth platform, the missile may be said to have arrived at its 
destination. 

The two pendulous units are mounted upon the controlled member, one 
with its plane of motion in that of the programmed great circle, the other 
with its plane of motion perpendicular to that great circle plane. The lon- 
gitudinal pendulum, i.e. the unit with its input plane parallel to the plane 
of the programmed track, sensing longitudinal accelerations, produces an 
output signal which is then modified through the implementation of the lon- 
gitudinal mechanization equation. The resultant signal is used as the input 





* This orientation is accomplished by rotating the controlled member on 
the earth platform to the correct meridian angle of the great circle track 
at the point of departure, [Ay (G.C.) (dep) > and by tilting the outer gimbal 


axis to the correct polar angle (PA) (dep) of the great circle track plane 
with respect to the polar axis of the ‘°°?’ earth. 
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to a servodrive system that rotates the controlled member so that it 
tracks the local vertical. Similarly, the track control pendulum, sensing 
lateral accelerations; (i.e. accelerations perpendicular to the great circle 
track plane), produces an output which is modified through the implemen- 
tation of the lateral mechanization equation, and the resultant is used to 
position the missile control surfaces so that the missile remains close to 
the programmed great circle track. The altitude of the missile above the 
surface of the earth is maintained by a pressure type altimeter that sup- 
plies the essential input to the altitude control system. 

A Coriolis acceleration computer is used to compensate for the ef- 
fects of the Coriolis acceleration upon the action of the system. Geodesic 
acceleration is in general small enough to be ignored. It will be the pur- 
pose of this thesis to discover whether it is theoretically possible to main- 
tain the errors between the indicated and actual local verticals small 
enough to permit sufficiently accurate indication of the local vertical for 
the long-range bombing and missile guidance problems. Various mechan- 
ization equations will be examined, to see which offer the greatest chances 
of success. The most promising of these will be furthur examined to de- 


duce the optimum control parameters. 
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CHAPTER 


KINEMATICS 


1, Purpose 

It is the purpose of this chapter to derive the fundamental kinematic 
equations describing the motions of the missile postulated in Chapter I. The 
principle geometrical relationships are illustrated in Fig II-1, and are shown 
symbolically in Fig I-2. 

The following simplifying assumptions are made in deriving the kine~ 
matic equations: 


Lateral and longitudinal motions decoupled 
Constant altitude flight path 
Spherical and homogeneous earth 


Controlled member located at center of gravity 
of missile 


e. Zero angie of attack and side slip angle 


a 0 Ff & 


The decoupling of the lateral and longitudinal motions divides the prob- 
lem into one of control of the missile in the direction perpendicular to the 
programmed great circle and one of indication of the distance travelled in 
the programmed great circle. The kinematic equations are therefore de- 
rived separately for lateral and longitudinal motions. 

Since the flight path is to be considered as having constant altitude above 
the surface of the earth, the only vertical component of accelerationis the in- 
ertia reaction acceleration of gravity. With the earth assumed tobe spherical 
and homogeneous, this componentis a constantover the surface of the earth, 
with the numerical value of 115,920 feet per minute per minute (equivalent to 
32.2 ft/sec”). The assumption of aspherical earth also eliminates any geo- 


desic acceleration terms from the kinematic equations (ref par 1 App A). 
an 
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A Cortolis acceleration computer is included in the fundamental guidance 
- system (Fig I-5), and is assumed to work ideally (ref par 2 App A ). The 
kinematic equations, therefore, do not contain the Coriolis acceleration. 

With the controlled member located at the center of gravity of the mis- 
sile, the gimbal isolation system can be considered to remove entirely all 
effects of the roll and pitch of the missile. No accelerations reach the con- 
trolled member as a result of angular accelerations of the missile about its 
body axes. 

Considering both the angle of attack and the angle of side slip as zero, 
the velocity vector of the missile always lies along the longitudinal body 
axis of the missile; that is, Vi (air) -M] lies along Xn (see Fig II-3). This 
is taken as the control direction (CD). The kinematic equations for lateral 


motion are derived first. 


2. Kinematic Equations of Motion Perpendicular to Great Circle Track* 
The geometrical elements of the simplified lateral guidance problem are 


illustrated in Fig II-4. Figure II-5 defines the directions and angles that are 
used in the derivation of the kinematic equations. 


As can be séen from Fig II-5, the following relations are true: 


(CVA vertXte) = (D Ch eruey(te) ~ (D CO) indy te) (11-1) 

Gre-(emy (te) z Re (CA (vert (te) (1-2) 
~ Die ccmilt 

tan(DC) ¢ruey(tcy = a (11-3) 


*This development follows that of John Hutzenlaub, Notebook, July, 1947. 
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x _~% ; > = ; 
foe Pp ie 
j Key f % 


j ia 


Mie-(cm)) ey) ViE-(airy) (ce ee m7 SinAz ¢ rm) cos(C) A Werty(tc) 


(II-4) 
From (II-2) and (II-3) 
ROA wertyitey = 7 919 tan (DOV it pue rte 
(II-5) 
But, from (II-1) 
tan(DC), veyto = tan VA Wert yiter . (DO) indytey! 
_ tan(O)A wert ite) + tan (DC) (ind vite) 
1- tan C)Awyertyitey tan(DO)indyitcy 
(1-6) 
Therefore, 
ma. _—_itan(O)AWvertyitey 
_ One (C) Awwert)(tc) 1-tan(C)Awyerty(te) tan (DC) cindytey 
_tan(DO)indxtey 
1- tan (C)A Werty(te) tan(DC) (ind)itc) 
(II-7) 


This is a theoretically correct expression, if no vertical accelerations are 
present, or if vertical accelerations are assumed to cause a variation in 
Erp: However, if the size of (C)A wert) (tc) does not become greater than 10 
milliradians, an error of no more than one part in 30,000 is introduced when 


(C)A is substituted for tan(C)A (vert)(tc): This should be true for all 


(vert) 
cases of practical interest. If this is done, eq (II-7) becomes: 
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. a. (C)A Wertxt 
CA... is eee SVEROXIC) Ol 
(Vert)(tc) ( Re 1-(C)A wertyite) tan (DC) (und) tcy 





(28) tan (DC) indyite) 
Re? 1-(C)Awvertyaey tan(DC) indyte) 
(II-8) 
Err E 2 
If Re is defined as W_,_, then eq (11-8) becomes: 
ae (COA t t 
(C : i Woe (Ver d( Cc) — 
Awvertxte) "™ 1-(ClAwertyitey tan(DC) indy tes 
_we tan (DC) (ind)(te) 
a (C) A werty te) tan(DO)indy(te) 
(II-9) 


If (DC) (ind) (te) is smaller than fifteen degrees,* (DC) (ind) (te) can be 
substituted for tan (DC) ( nd) (tc) for first order theoretical solutions. Also, 
(C)A (vert) (tc) tan (DC) (ind) (te) , as the product of two small angles, becomes 
very much less than unity, and can be neglected. Then: 


ee 2 2 
(CVA ertyitey * Wane C)Awertyctey = 7 We OO) inate) (II-10) 


Using the same small angle assumptions, and the additional assumption 
that Az (em) remains less than fifteen degrees, eq (II74) can be written in 
the approximate form: 


Vie-com)) (tc) * Vie-(airy) (tc) — V itair)- M] Axicm) 
| (1-11) 


* Ref par 3 AppA 
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Equations (II-2), (11-10) and (II-11) are the basic kinematic equations that 


are used in all the derivations of Chapter II. 


3, Longitudinal Motion Kinematic Equations* 


The longitudinal problem is one of range indication. The geometrical 


elements of the simplified range indication problem are illustrated in Fig 


II-6. Figure II-7 defines the directions and angles that are used in the der- 


ivation of the kinematic equations. 


As can be seen from Fig II-7, the following relations are true: 


(DC) (trueydongy = (C)A(vertylong) * (DC) cindy long) 


A = Qe (emy (long) 
r(true) R_ | 
Gi e-(cr) (lon ) 
tan(DC),4-ye(1ong#—— —— 


Gir 


Taking the tangent of eq (II-12) 


tan(ChA + tan(0C),. 
Rap (DO) t ruey(1ong)= eee ee none). 


1-tan (C)A (vert )(long) tan(DC) (ind) ( loaq) 


(I-12) 


(11-13) 


(IT-14) 


(11-15) 


Equation (II-15) is a theoretically correct expression but, as in the 


lateral case, all practical conditions find (C)A (vert)(long) 


milliradians. Then, eq (II-15) simplifies essentially to: 


smaller than ten 


*This development follows that of John Hutzenlaub, Notebook, July, 1947. 
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a. ( , % , . : t } 
yw e i § j ¢ 


Aili 


tan(DO)/:-uercionay = _C) Awwertyciongy * £47 (DC indytongy 
J I-(CVA (Wert long) tan (DC) (indy, long) 


(11-16) 
From eqs (11-13) and (1-14), 
A (t Re A, (T 
tan(DC) | = ri true) — rue) 
(true) ( ong) Ore Wie (1-17) 
Substituting eq (11-17) into eq (11-16) 
Ar(truey  _ (C)A Wert) (long) * tan (DOI ind (long) 
NE 1- Q)A (erty (long) tan (DC) (indy (long) (11-18) 


From this there results, as the general longitudinal kinematic equation, 


2 (CYA (Wert) (lon ) 


= a a Ne 
Av(true)= 


mE 1-CAwvert) (long) tan (DO) (ind) Cen 


A We Tan (DC) (indy(1ong) 


"A= (C) Aw erty (long) tan(DC) (indy (tong) (I-19) 


When the same small angle assumptions are made for the longitudinal 
kinematic relationship that were made in writing eq (I-10) in the lateral case, 
eq (11-19) reduces to 


A, (true) — Wie (C) A\vertylong) a Wne (BC) (indy( long) (II-20) 
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This equation and the relationships 


Ap(true) - Ay (ind) - (OVA (yerty(1ong) 


+ 


(1-21) 


and 


A = Ay (ind) + (DC) (indy(long) (I-22) 


r(app) 


form the basic kinematic equations that are used in all the derivations of 


Chapter IV. 
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CHAPTER III 
THE TRACK CONTROL PROBLEM 


1. Introduction 





The simplified kinematic equations developed in Chapter II 


for the track control problem are restated here for convenience. 


Ace (emt (cca Re (C) Avert) (tc) (O-1) 


(C) Awe rt)(tc) + W,, E ("= on —W.. (D C), ind)(tc) (OI-2) 


Vie ~(erm)] (tc) =Vee-caiv) (tc) —Vicairn-m3 Az cm) (1-3) 


Using these simplified kinematic equations, various mech- 
anization equations are investigated in an effort to determine a 
physically realizable and reasonably accurate mechanization of 
the track control problem. As previously stated, for the purposes 
of this study, the track control problem is considered as de- 
coupled from the longitudinal problem(1) | and vertical accelera- 
tions are neglected(2), It is further assumed that the missile 


operates with perfect control surface servomechanisms and that 





1 - ref par 2 App A 
2- ref par 4 App A 
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it possesses perfect aerodynamic response?) . Finally, the pendulous 
accelerometers are considered to indicate instantaneously the direction 
of the resultant acceleration”), 

The purpose of the track control mechanization is to reduce toward 
zero any angle (C) [ (Vert) (ind) ] (te) that has been introduced by lateral 
motion of the missile away from the prescribed great circle track, and 
to maintain this angle sufficiently small to insure arrival within one 
mile of the destination. The control is achieved through operations upon 
signals received from the measurement of the available angle 


(C) [(Vert)app |i ,4)> the pendulum angle, in order to correct the head- 


ing angle, “TX em) - (Cb) }° 


Throughout the derivations which follow, it is convenient to refer 


to Fig 11-5, which indicates the relationships among the important angles. 
2. Proportional Control 


The simplest possible mechanization is proportional control, in 
which the heading of the missile is changed in linear relation to the 
pendulum angle as measured between the pendulum and the controlled 
member. This case is therefore considered first. The mechanization 


equation for this concept can be expressed as: 


Aztec = — Steccoitaal (DC\indytey (1-4) 


(3) ref par 5 App A 
(4) ref par 6 App A 
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Using this mechanization equation, and the simplified kinematic 
equations, numbered one through three, a performance equation 


for this system is now derived. 


From eq (III-3), 


A == Vie inn (tc) ce-(cm)) (tc) 


2(cm) \V (11-5) 
C(/air)-M) 


Integrating eq (III-1), assuming the constant of integration as 


zero(1) | and substituting into eq (III-5), 


Vee. (aivyi(tc) Re E (CVA erty (EC) 


Ai« = (111-6) 
mM) 
A. MJ 
From eq (III-2), 
OA, ert)(tc) 
—(DC)inareo™ can HIC)A A wert yite) (I-77) 


W ve 


1 - par 7 App A 





Substituting (III-6) and (III-7) into the mechanization eq (III-&) 


Vee-ainae Re (Aware — 5 OAseee 4 (C)A | 
VV Ce(c DJ CAA) \/" Wert) te) | 
C(air)-M)J 


(111-8) 


From this, the performance equation can be written: 





R. Wu 


Vain -M] 1 [c(cO)J TAA] 


ae : | VV —_ 
, Ae a (CA evegt W.. (COA waa | Ee Ce-fair)j (td 


LAA Vicar. M] (1-9) 


» 


Although this equation provides a stable control system, a 
forced error results from the lateral component of the velocity 


of the wind, Vig 7 (air)) (¢¢)° This, obviously, is unsatisfactory 


for the control of guided missiles, because no large steady-state 


errors can be tolerated. 


3. Integral Control 


Next, let us examine a system that changes the heading of 
the missile in accordance with the integral of the pendulum angle. 


The equation for this mechanization can be written 


(141-10) 


z(cm) 


on c-h (DOV insets ot 


dae 





















Differentiating (III-6), and substituting with (III-7) into 
(III-10) differentiated, 


Vee-aienita Re OA versie 


= Fit DYTAA] 


(NA veces 


\A/- + (C lA Ve rt}i(te) 


Vea r)-M] 


2 ae 
(C)Avertyteg — 
}Fe + 


D. (co)J[AA] a . 


2: 
(C )A Vert)itc) a — 
g R + 2 Cc (co) LAA) V 


Ay 
ie 


This system, although the forced error is now due to the 
acceleration of the wind, has no damping term. S{(c)(cb)} [AA] 
can be made large to reduce the gize of the forced error, but 
this results in an undamped period of oscillation of approxi- 
mately 84 minutes for the controlled member. Any residual forced 


error (since it is physically impossible to make Sf(c)(cD)]} [AA] 


infinite), causes this oscillation to appear. 


4, Second Integral Control 





In order to reduce the forced error, consider next a system 
that changes the missile heading in obedience to the second inte- 


gral of the pendulum angle. Here, the mechanization equation is: 


ae [f Ch. (tc) dt dt 


Vv 


(}11-11) 


E- fair) (tc) 


(111-42) 








Differentiating (III-6) twice, and substituting with (III-7) 
into (III-13), differentiated twice, 


Vee —(air)j(tc) - Re. (C) ee r ccc ji 


> (C 
V tewi r)-M} 4 mere“) Wo. ) Ave rt (tc) 


on, 


*~ V, aiur)~M ‘, ae 
(C) A weces (tc) eS (c (CD)J [AA] A Pi iiecdiec om S Mca ae = V, 


Cc(c OL AA) es 
a. cOyLAA R. R. E+(AVratte) 


By Routh's stability criteria(l), this is unconditionally 


unstable. 


5. Derivative Control 


In order to give a complete presentation of possible types 


of mechanization, consider the heading angle to be controlled in 


s 


accordance with the derivative of the pendulum angle: 





eon can (DC) ina (tc) 


(141-14) 


(111-15) 


(111-16) 











Cor 


: 
4 
« 8 é kh a. il . 


Differentiating (III-7), and substituting, with (ILfI-6) into 
(III-156), differentiated, 


Ricca (tc) Re (C)A (vert) (tc) __ | C)A «ven +(C)A 
a [c(co)) [AA] \J (vert) (tc) C1 | ay) 


C(air)-M] NE 


From (III-17) 


a VV \W/ ~ 


{K(co)] LAA) vn [fair)- mj ee ontias 


CA seryegt W.. I+ 


VEIEtS) 


By Routh's stability criteria for a cubic, this is uncondi- 
tionally unstable. Moreover, there is a forced error due to the 


velocity of the wind. 


6. Proportional Plus Integral Control 


As shown in paragraph 2, proportional control produces a 
Solution which, although it has a forced error caused by wind 
velocity, is nevertheless stable, while integral control (para- 
graph 3) develops a solution that does not have any forced error 
caused by wind, but which also has no damping, and is therefore on 


nstability. Thus it might be expected that by com- 


1 equations, a stable system without forced error due to 















ng eq: (III-6) and (III-7) and substituting with 






| (III-7) into (III-19) differentiated, 






aa 
ce 
i. 


Vee -(air)” iGec) R. (C) A (Vert)(t0) __ (C) A\we rt)(+c) ; 
— Ss [c(co)J [AA ws . +( C) A 











(C)A Ver c 
sf ert) (tc) +CIA wert) (ee) 


Veein oo e (Verv)(te c(cp)} TAA) We (111-20) 


From eq (III-20), 





2 coy) [AA A 
a ae (C ) A (Vert) (tc) “1 


[c (cp) [AA] Wise: r)-M] W we 


W ie R. 


ge 
[pn\*p° | ; aa >; ] CAA) _ W 
(C)A werexteoF S BN ieee EW we cleo TAAL (COA (Vert) (te) = V ce = card (1tt-21) 


S te (e001 [AA] Seic ol [AA] Vecniry-m] 


By Routh's criteria of stability, this system is stable if: 


D3 CAA] and eas (cpjtaa) ore positive numbers, 
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ration of the wind. This suggests adding second: integral 





rstem just discussed. 
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cy (DC) ind) (tc) -- S fe (co)) [A (DC) naytes at dt | (li t-2a9 


i (III-7) differentiated twice, 
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C)A evar . (CA ' 
er c) (Vert) (tc) ( an 23) 
D (con CA 4) = +(C)A cence + Sieoneeh ae . T (C) A wered(ec) 


= (), Aw rt) (tc) + (C) A E 
= Dic teort& AY (Vert) (te) {| \W/ 
NE 


W we 


NE 









| ; 


| Devon [A'A} (C) > bert) oer 


r (c(co)] (A A] 


2 


W oe 
> Deemed wi CA NE Vee cin (111-24) 


(c(co)ILA A} 2 ) 3 
S Wwe (C Awerveyh S eae ce . VV 
(c(epyJ LA A] ([e(ccov] CK A) {c(eo)} CA A) Cfair)-m] 





| ~~ ic (cp) CA A ( A i 
Ss C Wert)(eey W 


~ (c(c oA Al 


eria for a quartic, this 





pitie have the same sign. 
Bitty scx all negative), 
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jt 
i. fs 
1 if 
{ i 
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, aw) FF ~~ | 
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he wind velocity enters into the forced 


Likely that a forced error of this nature coulc 
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large enough to be of practical importance. How- 







ea J > 
or should prove important, it is possible to go 
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- performance equations. One further case is con- _ 










1 Plus First, Second and Third Integral Control 
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fm. VC), nae) S ee (D C)., ‘Dlee) dt—> a ae ff DC),arue) At dt-9 [c(epuca’k Aly (D C), tnd) (4c) dt dt dt (111-25) | 


ng eqs (III-6)and (III-7) differentiated three 





differentiated twice, eq (III-7) differentiated 





into eq (III-25) differentiated three times, 


OA Vert) (te "A 
yr a (C) A entre 










| Bune 05 OAR. A (Cena 44 

C r (Vert)(t<) (Vert) (tc) 

- (<(cD))[A A) \a/ = Q)A eves er , oJ LA A] Ww? +(CAL Se +5, cleo} [aa] cat = C)A (Vert) (tc) 
Né NE 





(111-26) 












- NE 









Secon CAA] (C)A, : 7 < W ue V (| )-27) 
VertXé) 
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E-fir)) 


; (c(cOJ [A A] (c) Rn z Ww. = CA AJ 4 ee CA A} (C) A verso i W.. S A] 
e(co)}IA 


oo So , 
x [< (co){A A) fe (co)) LA A) [e(co)} LA AD 


(C)A,, { er Wy 
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ae ta A} > Ai Veen 








(C) A (Vert) om 
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nts stability criteria for a quintic, this 


uation is stable if: 


(a) All coefficients of the homogeneous 
equation are positive, 





(b) 
Decor ac \al 
S S NE 
[ (cd) [A A] WW a Ce(co)JC'A “A a | 
eo | FZ S \ 
[c(coMA AJ Ie + [c(co)]LA Al \ | 
Bc A) Vetwier a4 cena A) J : 
\ - 
Sea A) i ph Deeteontk 4) A) » at aici ce Al \W/ * Dc nets | 


S:. (c 0) 





Fi 4 Secon \- + ee | 
| | y NE NE 
as a. , ak (co) {A A] am (cDILA A] 
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NE 
One ae oad D)] C ‘A aM / Ce(eov ICA A} 


IR + 





: CAA] Vi. (air)-M) 
ie CAA) 


AA) eo. + Secconth '3 W. A ee a \ 5 CC(COVILA'K} 2 7 < 
[cicD) ILA A) “[lair)-mM] EC(eD. 


aie aa Decor [A A} WwW Deeceonts Sumo 3 Ye | =A g 


NE 
eugene Poy DILA A] 


7 S V . fp Mes . 
J "Cair)-m) fe(cD ILA ‘AJ [(air)-mM]J (c(co)ITA AJ celeoyTA AJ [c(CO)ILA A] (c(coyvITA AJ 


Cefe JCA A) 


‘These restrictions permit the sensitivities to be either 
_ or negative, provided that they all have the same sign. 
, is of interest to notice that in each stable performance 


1 |) ee 


the forced error has the same type of coefficient, with 





' of the wind derivative increasing from velocity for 
as the perform- 


der 


es contro! to any derivative desired, 


ion increases in complexity. This coefficient is of 
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The forced error is thus Seen to decrease as the missile airspeed 


increases, and as the Sensitivity -{C(CD)} [AA] increases. Theo- 


retically, if this Sensitivity increased to infinity, 


there would 
be no error, either transient or in the Steady state, caused by 


the action of the wind. This is similar to the conditions found 


with integral control, eq (III-12), where, as the Sensitivity 


> (C(CD)] [Aa] increases to infinity, both Steady state and trans- 


lent wind errors reduce to zero. 


9. Track Control Correction as a Function of Wind 


Since the mechanization eq (I11-22) for proportional plus 








first and second integral control appears to provide a stable 






System without excessive forced error, this system is considered 






further. The performance equation for this system, eq (III-24), 





can be expressed in operational form: 





Dp (111-28) 
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If (C)A vert) is expressed in minutes it is a measure, in 
nautical miles, of the transverse linear distence of the missile 
from the desired great circle track. The vresence of the opera- 
tional symbol in the numerator of the equation indicates that 
there is no error, in the steady state, from an accelerating 
wind. If it is desired to find the missile response, Az(em), to 


an accelerating wind, it can be calculeted if eq (III- 28) is 
AZ( em) 


(CYA vert ) ( tc) 


multiplied by an expression for The result 


yields the heading of the missile as a function of wind velocity. 
Manipulating eqs (IIJI-1), {III-2), (III-3), and (III-22), and 


expressing the result in operational form, 


ie CAA) R. 
A, (cm) a Y IR 
(C)A ct) (tc) 
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Combining eqs {I1]1-28) and (III-29) gives, 
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The negative sign shows that the heading AZ(cm) of wie air- 
plane is correct to oppose the tendency of VIE ~ (air)] to blow 
the missile off the great circle course. This is true even if 


the sensitivities are negative. 


10. Track Control Closed Loop System 
The relationship of eq (III-30) can be obtained from a study 


of the over-all lateral system as a servomechanism. In the direct 
servo closed-loop system, the angle (C)A( vert) does not appear, 
but it can easily be obtained, as indicated in the accompanying 
Fig III-1l. 


In this figure: 


(PF), has the value b 
BIR 





(PF )> ia moe and is made up of the sum 

(DC (ind)(tc) of two performance functions, 
(PF). . (DC) ind) (tec and . 

DC) (true) (te) 


a , 
(PF)p = —Z(cm) Ss which is the mechanization 
(DC) (4nd) (te) equation in operational form 





(PF)3, in the feedback path, is (- ~lledr) = Mp, 


providing the proper function for the com- 
binor, which is physically the pendulous unit. 
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The pendulous unit, if (DC) (true) remains smaller than about fif- 


teen degrees, solves the equation 


wc A _e -(air)) (Dp C). 
g | Z(cm) a. (+rwe)(tc) ( i 3p 


This equation is derived from eqs (III-1), (IIT-2), (III-3) 
differentiated, and the relationship that 


(DO) cenucted= (C)A (Vert (te) 7 (DO), aye Ib b=32F 


This last equation can readily be obtained from an inspection of 
Fig II-l. 

Solving for the performance function (PF),, using eqs (III-1), 
*) (III-3), and (III-32), 


(DO)naree _ p +W\, 


(1 1339 
© (DC), .se9 rue) (tc) p : 


(PF),= 


The mechanization equation, which in this example is eq (III-22) 


expressed in operational form, is: 
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The product of eqs (TII-33) and (III-34) gives 
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but, from the theory of servomechanisns, 1) 


Nm _ (PF) (PF) 















ae | ++ (PF) (PF) (111-36) 
Solving this equation gives eq (I1I-30), as before. The stability 
criteria for eq (III-30) are identical with those for eq (III-2h). 
Since the sensitivities can have either sign, a necessary condi- 
° EIR 
tion for stability is that | Sfcicp laa | | 
if the sensitivities are chosen as negative. If this is expressed 
as 
QQ — ** : 


fe - ir 









"Principles of Servomechanisms". 
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then, for stability, 


-— sO | (111-38) 
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11. Numerical Solutions, Using Two Quadratics 

Using the relationships that have been established, and 
substituting numerical values for the coefficients, it is possi- 
ble to solve for the various angles as a function of the accelera- 
tion of the wind Vie - (air)]° From a sequence of such solutions, 
using a range of values for the sensitivities, it is possible to 
select the most satisfactory values, and to estimate the effec- 
tiveness of the system. This is now done. 

The homogeneous portion of the performance equation, (III-24), 
is a biquadratic. It can then be written 


(p'+ 2(DR\W, +) (p+2(0R), W, +) = 0 (111-39) 


which factors the equation into two quadratics. Manipulation of 
the equation indicates that if the damping ratios and natural 
frequencies are varied, the higher damping term associates itself 
with the shorter of the two frequencies. This result then appears 
essentially as that of a lightly damped long period quadratic, 
which is undesirable. For this reason (DR), and (DR)o , and W) 


and Wo have been chosen equal for all solutions of the equation, 
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and equations are now derived to facilitate the selection of 
values for the sensitivities to make these limitations hold. 


Expanding (III-39), in the general form, gives 


p*+[2(0R)w,+2(DR\ wi) p*+ [w.'+W, + 4DRYDR) Wi W,]p* + 








[2(DR) W'W,+2(DR)W,W,]p.+W? Ws, =O _ 
Wy Wo 
If tne ® We = (FR), and (DR) = (DR)> = (DR) then, 


comparing the coefficients of eq (III-40) with those of eq (III-2k4) 


gives the following relationships: 
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When these four equations are solved simultaneously for the 


sensitivities and for the damping ratio, the results are: 


(DR)= = | (111-45) 


Fg 


Sreteontk I= BW. Vin... (DRYER) (111-46) 
— 9, (FR) 
= A) 2Vecerag (DR) (111-47) 


5 = Fa Wa (FRY 
[c(e oy] CAA] BV 4.9. _ (DRY (111-48) 


It is noticed that the choice of a damping ratio immediately 
fixes all sensitivities and the frequency ratio (or that the choice 
of a a kiekey ratio, alternatively, fixes the damping ratio). As 
the frequency ratio increases; i.e., as the quadratics have pro- 
gressively higher natural frequencies, the damping ratio simul- 
‘taneously increases without limit. Conversely, as the frequency 
of the quadratics approaches that of the 8& minute period "earth 
pendulum", Wyp , the damping reduces to zero. The sensitivities, 


at the same time, increase without limit. 
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l2. Numerical Solution Using Four First Order Terms 
A second method of writing the quartic equation, which can 


be examined easily, and which may yield results of practical 
interest, is that in which the equation is reduced to four decay- 


ing exponentials. This equation is of the form ‘ 


PreallP en) ea) Pren)=° om 


When expanded, this becomes ; | | 
ae ere I ee et LI -—— 
a en) ten? ent er) | P+ lew on), Ver), ee,” ,* He, 





1 2. | | 
“, cay, |? P Ravan * (r) (cr) (er), * (CT) (CT), (cr), (ct), (ct), ( en |P 
“(on (en, Kor) (cu), O : (111-50) 
The general relationships among the coefficients of eq (III-50) 
and eq (III-24) give 
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The special case which is considered here is that for which 
(CT), = (CT) = (CT)3 = (CT), = (CT). Under these special circun- 
stances, the simultaneous solution of eqs (III-51, 52, 53, and 54) 


yields 
ul =e ( 
Cry (c+ +v32), (111-55) 
from which 
a 
; O. : 
(CT) =W.. ZAl8B o a =W., 412. (141-56) 


g Wwe 











Cele nyICA x Vices = 4W. ‘ 4 Chl 1-57) 
(CT) (CT) 
and 
ce > (CT) (111-58). 
~ (coLAA) >. a 
CQir-m] 2 





“« CT) 


a i 





CONFIDENTIAL 


S _ he (CT) 
» A re(eo)} CAA) VV wk 


Clair)-mj if 2 
7. =. (111-59) 
“ (cTy| 





Evaluating (III-57), (III-58), and (III-59), using eq (III-56) 
with the smaller (CT), 
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This gives negative values to the sensitivities. The alterna- 
tive evaluation would make the sensitivities positive, but would 
result in characteristic times far longer than could be tolerated. 
The system, with negative sensitivities, is stable by Routh's 
criteria. 

13. Numerical Solution Using Quadratic and Two First Order Terms 

The final method of examining a quartic equation which may be 
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of interest is that in which the quartic is broken up into one 


quadratic term, and two decaying exponentials, 


(p'+2(DR)Wp + Ww’) (p+) (p+ aa ll 


Upon expanding eq (III-63) this becomes, 
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The general relationships among the coefficients of equations 


(III-64) and (III-24) are 
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These equations can be solved simultaneously for the sen- 
sitivities, and for the relation which connects (CT), , (CT)o2 
ae (DR). It is convenient to use the relationships: (FR) = 


CT)4W CT) 5W 
( NE (CT) 2Wne at 


and {(CT)PR] 2 — 


= {(cT)PR} = 


esulting equations are: 
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These are evaluated for the single case of (DR) = iz » with 
((CT)PR], = [(CT)PR]> = 0.1. Then, from (III-69), (FR) = 0.307. 
This gives: 


5 7 Di. 
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Notice that once again the sensitivities are negative. 
14. Plots of Results 
It is now possible, using various coefficients for the 
equations of paragraph ll, or the indicated coefficients of the 
equations of paragraphs 1l2 and 13, to obtain numerical answers 
by solving the equations with various wind inputs. This was done 
on the Rockefeller Analyzer at Massachusetts Institute of Technology. 
It is not necessary to solve for various missile speeds, since 
the speed of the missile can be combined with the heading angle of 
the missile AZ (cm)° This is possible because of the linearisation 


assumption which was made in deriving the kinematic equations, and 
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holds only if the missile airspeed is very much greater ‘than the 
velocity of the wind. 

Plots were made of the angular error in minutes of arc of 
the position of the missile perpendicular to the great circle 
track. Considering the earth to be spherical, this is directly 
@ measure of distance from the track in nautical miles. For 
convenience, on all plots, the unit of angle is the minute of 
arc, and of time, the minute of time. The foot has been taken 
as the unit of distance, 

The first input of wind which was used in solving the equa- 
tions was a rectangular pulse of wind acceleration, or a constant 
acceleration of the velocity of the wind from zero to some finite 
value. The duration of the pulse was taken as one minute. Next, 
responses were taken for a step function input of wind accelera- 
tion. Finally, responses were determined for sinusoidal inputs of 
wind velocity using several frequencies. Plots of the results 


follow, using various sensitivities in eq (III-22);: 
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NUMERICAL CONSTANTS USED IN OBTAINING 
RESPONSE OF TRACK CONTROL SYSTEM TO PULSE AND 
STEP OF WIND ACCELERATION AND CIRCULATING WIND 
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CHAPTER IV 
THE RANGE INDICATION PROBLEM 


1. Introduction 
The simplified kinematic equations developed in Chapter II 


for the longitudinal problem are restated here for convenience: 


(DC) indy(onei™ —— (C) A\ Wert (long) 


A 
W ie 


(1V=1) 


Avcina)™ Av terue) + (C) A r (1V-2) 


) By tread + (DC). na) ions) (1V-3) 


Using these simplified kinematic equations, various mech- 
anization equations are investigated in this chapter in an effort 
to determine a physically realizable and reasonably accurate 
mechanization of the range indication problem. The assumptions 
are identical with those for the lateral problem. The missile is 
assumed to move with instantaneous response to wind forcing 
functions. 

The purpose of the longitudinal mechanization is to reduce 
toward zero any angle (C) [(Vert) (ind) (long)} which has been in- 
troduced by motion of the missile along the prescribed great 
circle keer, and by measuring the angle between the reference 


vertical and the indicated vertical, to obtain a measure of the 
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the missile. The orientation of the controlled member 
echieved by measuring the angle (DC) (ind) (long) between the 
9sition of the pendulum and the controlled member, and by oper- 
iting on this quantity to position the controlled member with 
respect to the reference vertical. The longitudinal reference is 
shosen as the vertical of the departure point, (Vert)(qep)- Since 
she value of (DC)(ind)(long) 18 dependent upon the position of 
the controlled member, as well as of the pendulum, it is more 
sonvenient, in deriving mechanizations, to use the angle Ar(app) 
between the reference vertical and the indicated vertical of the 
ontrolled member, than to use (DC)(ind)(long) directly. Ar(app) 
ig physically measurable, since it is the sum of (DC)(jnd)(long) 
and Ar(ind), Which are both known. It has tne advantage of being 
independent of the position of the controlled member. 
Throughout the derivations that follow, it is convenient to 
refer to Figure Iil-6, which indicates the relationships among 


he important angles. 


Trial Mechanizations 
A. The first mechanization that is considered is that in 


second integral of Ar(app) moves the controlled member 


cording to Ay(ind) and its first two integrals. That any simple: 


= 


oduces an unsatisfactory solution is indicated by the mec! 


equations for the track cortrol problem. The mechanizeé 


es /wp dt +S... = |] Aeon dtdt 
|“ cicemy tAa)  / A (ind) [c(em’ ha r (App) (1V-4) , 
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€ Fs iu 
ay - - 








Differentiating twice, 


-_™ Peas a [AA] Av cinad r a [AA] A vind) A vcapp) ( v=) 


Substituting (IV-3) into (IV-4) 


Seimatns | Ay inay Piers.) Ayeanay Drelemdti A) A indy = DC iinadiong b 
( 1V-6) 


which is the mechanization equation in terms of (DC) (ind) (long) 


and An(ing)- Further substituting (IV-1) and (IV-2) into (IV-6), 


| | A Pe. Pia A _—_ i as » | A =(C)A + — C)A, a Sreennen (C)A. 
mn! r (true) r (true) S S r{true)_ r S S 
W Die MILA Al conti A) [c(Cm)ILA A) [c(cm)} LAA} [c(cm)] [A A] [c(cm)} [AA] 
(1V-7) 
This control equation indicates that the value for S{c(cmM)] [aa] 
should be unity, and it will be so chosen in further considerations 
of the problem. Even with the position error eliminated by the 
proper selection of S{(c)(cm)] [aa] , there remains an error caused 
by the rate of change of Ay(true)- Since =ritrue) i3 a measure 
E 
of the distance travelled over the surface of the earth, and because 
airspeed is constant, 
A vere = 
p oo 
€ ( 1V=3) 
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This means that the steady state error is a function of 

wind velocity, and that this equation is therefore unsatisfactory. 
B. In order to improve this response, the first integral of 

Ar(app) 18 added to the second integral in the preceding mech- 


anization equation. The resulting equation is 


Jp ring At dt + Dreier th 3) dt +S {(cm)} TAA] 4) a Ay app dt dt tr ae A aoe) dt Ve) 


from which, with (IV-3), 


emai. Aviind) a (DC). sion a (DC). song) 








CIVH-10) 
Further substituting (IV-1) and (IV-2) into (IV-10), 
| ie riarka, . a, | — | | A =z (CA i Deeremnacaad (C)A A 8 (C)A. 
A vere r(true) 7 ‘ ( I\V-11 ) 
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In this control equation, there is a forced error resulting 


from the acceleration of the wind which can be removed if 


is made equal to Wwe. Under this condition, a 


trol equation results with a forced error resulting 
t of the wind of lower order than the time rate of 
sration. The homogeneous equation contains damping, 


4 minutes, because the natural fre- 
==: 
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as | Wee wh the value of S [o(¢m) [AA 























Because it is desirable to increase tris natlral free 


nN order to reduce the time during which transient errors 





rtance, otner mechanization eguations are examined in 
feet 2M an effort to acquire control over period as well 
ing in the control equation. 
The third mechanization equation to be considered uses 


d integral Ot Ar (app) to move the controlled member, by 


ating on Ar( ind ) and its first tnree integrals. This mech- 


mn can be expressed 


[//A, 4 dt dtdt + ee ae Ayeiine dt at+ Ser m) J [AAT Aye ing) at q Sec: <m)]C AAJ A,, nd) Lp "(9 PP) at d 


(1-12) 


rentiating equation (IV-12) three times, and substitut- 
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tion (IV-3) into (IV-12), 
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fc) -13) 
substituting ezustions (IV-1) and (IV-2) into eq 
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(1V-14) 
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In an effort to improve this situation, the second in- 


tegral of Ar (app) 1s added to the third integral in the preced- 


ing mechanization. Then, 





' Tf bse dt dt dt + —— aw | [A r(ind) dt dt a CK AY ind) BETO ccccmarite Arf Beare dt dt dt i: ee 


from which, with eq (IV-3), 


Bescon Pr tindyt Stefernt CAA) A tind) = (DC);i.4)10ng " Piecemucie (DC 


Further substituting eqs (IV-1) and (IV-2) into eq (IV-16), 
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The lowest order of the forced error in this control equation. 


is that of the acceleration of the wind. The coefficient of this 


S e@ 
erm, however is | (c(cM)] f AA) ne |: ner 


S(C(CM)] {AA} = Wee S(c(cM)]E AA] 


C(CM)} [AA] is chosen age equal to - , the acceleration term 
- W“NE 


ppears and the control equation indicates that the lowest 
| forced error is that resulting from the time rate of 


the acceleration. The coefficient of this term is of 


to that of the acceleration term , and is 





1: CAA} (C)A shes 


De (CmyyCAA] 
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Le(Cm]L AAI 
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Cc(cm)] CA‘ 


Axe at dt 


(1V-15) 
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ve S{c(cm)] [AA] 


. It thus appears that the acceleration 


5 
erm can be made to disappear if ic cM L AAJ is chosen 
[C(CMYJC AAJ 


NE ° 
From Routh's stability criteria for a cubic, the sign of 


term in the homogeneous equation must be sositive, and 


cmnyI CAA) tecomtiay 
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mYICK AY [c(cmy CX A} S [C(CMIILA A) 


tion (IV-18) shows that, with Sfo(om)][ AA] = — » the system 
NE 


table only if “tc (cm)] CAAI 


— is greater that Wy, . The 
> fc(cm)] [4a] 


tleration rate term can therefore not be eliminated if the 
mis to remain stable. The equation that has been considered 
_if it is physically possible to mechanize it, appears to 

de a satisfactory system. It remains to be seen whether 
[onal terms in the mechanization ecuation can improve the 

E. Consider next, therefore, the performance which results 


-o the mechanization enguation just considered, the first 





‘al of Ar(app) is added to the second and third. This mech- 


| nd) dt dt dt a on (cm)ICAal | A r(ind) dt dt .. m)} LA AJ A CH 


a 
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nich, with eq (IV-3), 


















8 (1V=20) 
tc (em) (KA A cindy = (DC); (ind) long + Si (cm) LAA} (DC), ind) long TS Cet metind) tong 
wther substituting eqs (IV-1) and (IV-2) into eq (IV-20), 

a h = a... — StecenmIL hay A +- aa | A =(C)A i Preientin CYA + > recom cial (C)A + i (C 
12 A r (true) ZZ \W? r (true) Ven r (true) r _ r S : r 5 VA r (1V-21 ) 
Ne Dyess [AAJ NE Jretemy [Kal = emyJLA AJ [c(cm)JLA Al [c(cmiILA A] (c(Cm)] LAA] 





re again, the lowest order forced-error term contains the 


The coefficient of this term is now 





ration of the wind. 


and does not contain any controllable sensitivity. This 


is therefore less satisfactory than that of par 2-D. 








For the sake of completeness, Ar(app) is next added to 


grals in the mechanization equation just discussed. This 


MICA Al A x cind) = H/o dt dt dt Tee ential coer dt dt =... ree A... pp) at ? ae CA A) Aap) 
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Further substituting eqs -(IV-1) and (IV-2) into eq (IV-23) 
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r (true) A. ruc) 
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$ system, like that of mr 2-E, contains an unremovable 


ly~state error proportional to the magnitude of the wind 
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When examining possible mechanization equations, it 


seems reasonable to add to the system of par 2-D the fourth 


~4 ue : _ 
gral of Apn(ing), using the fourth integral of Ay(app)- This 


NE S (c(cmy]T A A] 


A iru) = 


Ae r(ind) dt dt + De. (cmy)L AA} A. ind ) at + S [c(cm)J TA Al A, nd) =|f//1. G : pp? at at at dt 
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As was true with the corresponding system of par 2-C, an 
or caused by the velocity of the wind results. 
H. When the third integral of Ar(app) is added to the 


rth integral used in par 2-G, the equation becomes 


PfArins dt dt dt dt + Semin [/[Arinn deat dt + 


m which, with eq (IV-3) 


(cmLAA] Avtindy = LAAI A, nd) Sree’ A. ind) (DC), ind) long r 


5 S 


= Ec(omJCA al _ Cc (Cm) CA AI A ae 2 Cc(cmy} CAA) 
rue) 
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Z 5 Wye 5 Cele (A AY S 


(c(cmyJLA A) CC(Cm)] LK al 


Here, as with the system of par 2-D, the control equation 


“ates that the lowest order of error in the steady state 





s from the acceleration of the wind. The coefficient of 


S{c(cm)] (4a) Wenge S{c(cm)) (Aal 


C Aa] p A Seeiem RA A (core k ‘ mJ CAA) 7 
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(1V=23) 
cc(cm)JLAA (DC) (ind ) long 
(1V-29) 
‘A A > Maiti A Seine / 5h. WIAA] , | 
A r (true) (C)A, - (C)A,+ 5 : (C)A+ . (C)A. i YS (C)A, 
cm)}t AA] eee [c(cmi]L A A) Sectemit¥A [ce 1A Al (1V-30) 
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As before, if S [c(cm)} [AA] is made equal to ae the 


acceleration term disappears. The coefficient of the acceleration 


Stc(cm)]£aa}] — S{c(cm)}Caal 
Stc(cm)) (Aa? W*ne Sec(cm)]C aad 


rate term is It appears 


S 606 
that this term can be made to disappear if ?tc(cM)] CAAT is 
S(c(cm)} [ aa J 


chosen equal to = ‘ 
W°NE 


From Routh's stability criteria for a quartic, the sign of 


each term in the homogeneous equation must be positive, and 


2 2 
= are Seen [AA] Sete CRA Mee itie }~{S.cauera > 0 (1V=-31) 





If STc(cM)] [AA] is made equal to ~ , eq (IV-31) becomes 
W°NE . 
S S 7 Z C Z 2 
[c(cm)} CA A] Cce(cm)] CAAT | - [c(cm)C AA] 
Dectem tk es Daerah > O 


Bens CAA] W 


S cctem ( IV=-32) 


Therefore, in order to eliminate the acceleration rate tern, 
= Stc(cp)}{ AA] [coe] Chal} * must be greater than zero, which 
is possible only if S(c(cpj}[ aa] becomes negative. Routh's first 
Criterion of stability, which forbids variation in the signs of the 


quartic coefficient, is then violated, and the system is unstable. 


at 


















t has not been possible to eliminate, 
effect of wind of higher order than 
ition eq (IV-28) appears to give the 
steady state, that results from the use 
sing between these equations, it is 


response and ease of mechanization. 


et to. el iminate the effect of acceleration 








: 1 of Ar(app) i$ next added to the mech- 
par 2-G. This gives 


eT Bee in | A r(ind) at at at Tee (cm) A 


SFr) (DC),,narieng tS ce(em) JC AA) (DC),,né)10%g 
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(IV-1) and (IV-2) into eq (IV-34) 
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e the acceleration rate term is 
= a—£ 


nte aol *, avn 
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) and of Ar(ind) to 







in par 2-G. This gives, 





off r(ind) at dt at + Sema fA -tina) At at +S cent ins at T ee Aen fff r(app) at dt at at at Scansialff A r(app) dt dt dt at 


( 1V~-36) 
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J-1) and (IV-2) into eq (IV-37), 
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th's stability criteria for a quintic, the co- 


of the homogeneous equation must all be positive, 
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The system is stable, with proper choice of the other co- 
efficientsyif Stc(cmyl{ AA] is made equal to a eliminating 
NE 
he acceleration term. If this substitution is made in eq (IV-39) 
S eee 
nd if in addition the substitution ic (cM)It AA] — is 
SIC(CM)][ AA] WNE 
2. 
S _. aia, CA a) 
“Tclemyits GS a 
— si J celemaUKay a - | NE Ww. e 
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Since all of the sensitivities are positive, this equation 


2 
(2) -3 m 
can be divided by we ic(CM)JCAA) | arovided that 


Stc(cm)]{ AA] 


ee O¢ 2 
only S[c(cD ai AA] is maintained smaller than ( 52 : 


Without changing the inequality sign. This yields 


the | z 








Drctemyricn cin ii A Al \a/’ - Deve myIC AA \n/* 


AJ A 


Or, simplifying, 


. _= CA AY > O 


Since Spccm)}{ aa] must be kept positive in order to satisfy 
the first Routh criterion, both criteria cannot simultaneously 
be satisfied, and the system becomes unstable if an attempt is 


made to eliminate the acceleration rate steady-state error. 
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3. Implementation of Equations 

A. Four mechanization equations have now been discovered 
that eliminate steady-state error caused by the acceleration 
of the wind. It has not been possible to eliminate any higher 
order terms, but is seems unlikely that any steady-state error 
resulting from the acceleration rate of the wind can have appre- 
clable effect. The simplest of these equations was eliminated 
because natural freouency could not be controlled. The three 
usable equations are (IV-15), (IV-28), Bad (IV-3€). In eq 
(IV-15), two of the three sensitivities can be controlled; in 
eq (IV-28), three of four sensitivities can be controlled; and 
in eq (IV-36), four of five sensitivities can be controlled. In 
order to obtain the most satisfactory possible dynamic response, 
the use of one of the more com-licated equations may be indicated, 
but this choice of equations is subordinate to the considerations 
of ease and accuracy of mecnanizétion. A practicable method of 
implementing eqs (IV-15) and (IV-28) is now presented, and the 
complexities compared. 

B. In mechanizing eq (IV-15), (DC)(ind)(1long) is used 
instead of Ar(app) as the output of the pendulous accelerometer, 
as explained in par 1 of this chapter. The resulting eq (IV-16) 


fs restated here for convenience, in integrated form: 
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Now, an imaginary auxiliary direction is established, as 
shown in Fig IV-1. This direction is controlled, with respect 


to the reference vertical, by the equation 


PY nurvret= SD cceaunataiy [2.x long dt at 


(1V=45) 


The controlled member is positioned from the reference vertical 


through the use of the auxiliary direction by the equation 


A indy — Drecem)auxd tak: [00s at a sy | Arciaddiate dt (1V-46) 


These equations, when combined with the relation 
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This is identical with the desired eq (IV-44), with the 
relationships among sensitivities that 
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AUXILIARY DIRECTION FOR CUBIC MECHANIZATION 
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Equations (IV-45) and (IV-46) can be mechanized easily, as 





hown in Fig IV-2. The mechanization of eq (IV-15) is therefore 
easible. 
C. Equation (IV-28) is now to be mechanized in a manner 

milar to that used for eq (IV-15). The integrated form of the 


yuation, using (DC) (ind) (long) instead of Arjapp), as before 
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to implement this equation. 
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general closed-loop equations are derived in this section. 


The general mechanization equation can be written, for 


this study, as 


A vind) = (PF Lae ee long |V-66) 


The linear distance of the missile measured along the track, 


from the initial reference point, can be written as 


X= []Nec-ain dt dt TN cea aoa + X dinitial ( \V-67) 


Also, 


Ar =p 
P | r(true) (1V-68) 


E 


These relationships, with eqs (IV-1, 2 and 3) and the small- 


angle pendulum equation 


(DC) Bey. Aetendie Vee. (eins 
(true)lo ng °g a a} ( 1V-69) 
IR IR 


combine to form the closed-loop diagram of Fig IV-6. 
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The initial conditions in eq (IV-67) do not affect the 
response of (C)A; to Vig -~ (air)] » and can be taken as zero. 
Figure IV-6 can then be simplified to Fig IV-7. 

In accordance with servomechanism theory, the overall 
performance equation of the range indication system can be 


written 


(DF)... = (C)A, os (PE) sendy (PF neck) — (PP cind 





a I+ (PF), oy a 
where (PF) ( pend) = z and (PF) (kin) . — ae » rhe 


mechanization equations to be examined can be written 


EF Seecomncias Po (1V-71) 
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(PF) Speco P 


(mech) ( wartic = ; ; 
4 >: p + telemuten Dp + Jee commie p 
(1V-72) 


Using these equations, numerical values can be assigned to 
the sensitivities, eq (IV-69)solved for various types of wind 
forcing functions, and the responses using eqs (IV-72) and 


(IV-71) compared. Numerical values are first assigned to the 
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cubic mechanization eq (IV~%i), and then to the quartic eq inl 


5. —— of Numerical Values 


A. The homogeneous control equation for the cubic system can be fac- 


tored into a first order and a quadratic term, becoming 
(p+ —)(p*+2(OR)Wp+W')= 0, av-13) 


which expands into 
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Comparing coefficients term by term, with 
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Then, for any desired characteristic time and quadratic natural fre- 


quency, a damping ratio and set of sensitivities can be computed. Notice 
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that, once the characteristic time and the natural frequency have been 
selected, no control remains over the damping. | 

B. The homogeneous control equation for the quartic system can be fac- 
tored into two quadratic terms. As in the track control system, these are 


made equal, giving, 
(p +2 (DR)Wp +W’)(p+2(DR)Wp +W’)=O (IV-78) 
This expands into 
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Equating coefficients with those of the quartic control equation, with 
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Plots follow for responses of the cubic and quartic control equations, 


with various sensitivities. 
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CHAPTER V 
CONCLUSIONS AND SUGGESTIONS 


1. Disturbances to the Control Systems 


In examining the success of the mechanization systems whose responses 
have been analyzed in the preceding chapters, it is necessary to consider the 
types of disturbing inputs which will affect the systems. These disturbances, 
or winds, are of three major types. First, there are gust disturbances of 
fairly high frequency and short duration. These are enorma@usly attenuated 
by the long time constants of the systems analyzed, and should cause no 
trouble. Second, there are acceleration pulse type inputs, such as would be 
encountered when crossing weather fronts, where the wind velocity changes 
very rapidly to a finite and fairly constant value. These can be analyzed by 
the step function responses of Chapters III and IV. In Fig V~-1, the path 
through the "eye" of a hurricane has been approximated with step functions 
of acceleration, and in Fig V~-2, the response of the track control system is 
shown to be satisfactory, even with this violent disturbance. 

There is, however, a third type of disturbance; that in which a wind of 
approximately constant velocity slowly rotates.* This condition has been 
approximated by sinusoidal inputs as shown in Chapters II] and IV. The error 
if the sinusoidal forcing function continues for a long period of time, is seen 
to be very large. A wind of about forty knots, rotating with a period of about 
one hour, will introduce, in the steady state, deviations from the track, in 
the track control system, of about 20 miles. 


*Ref. A. G. Bogosian Notebook dated March 18, 1947. 
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FIG.-Z%-2 RESPONSE OF TRACK CONTROL SYSTEM TO IDEALIZED HURRICANE. 
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Three ways of correcting this condition suggest themselves: 

a. Build systems with variable coefficients, which provide the neces- 
sary high damping with step-inputs, and have a different response to such 
low-acceleration inputs as sinusoidally varying winds. 

b. Increase the order of the track control system so that the forced 
error results only from high order greatly attenuated terms. 

c. Modify the system to reduce the effects of the forced errors. 


These three methods are all examined in the next section of this chapter. 


2. Suggestions for Improvement of System Response 
A. Track Control 


a. A variable coefficient system can be of value only if changing the 
damping and period of the system will improve the response. Since the co- 
efficient of the acceleration rate forced error term increases rapidly with 
increase in damping, it appears that this should be possible. A seriou of 
track control responses using the same magnitude of disturbance but dif- 
ferent damping ratios has therefore been examined, Fig. V-3. It is seen 
that no choice of damping is satisfactory. If damping is entirely removed, 
the forced error of the simplified system reduces to zero. Actually, a 
forced error will still exist if there is any vertical wind acceleration, or if 
non-ideal components are used. Then, if the frequency of the rotating wind 
is equal to that of the system, very large errors can result. This approach, 
therefore, does not provide the needed improvement of response. 

b. The purpose of this section is to investigate the response ofa 
quintic performance equation. From eq [11-26 the performance equation 


may be written: 
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PERIOD OF (C) A (vert) (tc) (max) 


IN MINUTES OF ARC 
(DR) ROTATING WIND [EQUAL TO NAUTICAL 


MILES OFF TRACK] 


1.875 10 71.735 
1.875 15 14.310 
1.500 15 13.620 
1.050 15 7.600 
1.875 20 13.033 
1.875 30 9.548 
1.500 30 14.899 
1.050 30 13.652 
1.500 40 12.080 
1.050 40 14.096 


MAXIMUM TRACK ERROR FOR VARIOUS (DR)s WITH ROTATING WINDS 
OF 31.4 MPH 


FIGURE V-3 
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study it is assumed that the denominator of the performance 


yuation may be factored into the form: 
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Then expanding and equating coefficients of like items, 
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Examination of eqs (V-2, 3, 4, 5, and 6) reveals that if K and (DR) are 
selected, then CT and all the sensitivities are fixed. General expressions 
are now developed for determining the numerical values of CT and of certain 
ratios of sensitivities. 


From simultaneous solution of eqs (V-3 and 5) 


4K (DR) (K*~ 1) 


Noel 2S on 
Wwe (-K* + 6K*- 1) (V-7) 


From eq (V-6) 
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From eqs (V-4, 6 and 7), 
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Limited time permits the evaluation of only one numerical solution of 


the performance equation. 


Let K = 2 andr 1. 
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(C) A (Vert) (tc) = .8 minutes of arc (.8 nautical miles) 

Comparing this answer to those obtained in the plot of results in Chap- 
ter III indicates that the quintic performance equation produces a better re- 
. sponse than the quartic one. It should be mentioned that no effort has been 

made to select optimum parameters for the quintic. 
| c. It is seen, by an examination of the track control system perform- 
ance function, that no control over the gain setting exists. As was shown in 
Chapter Ii, it is desirable to have the sensitivities negative for the quartic 
performance equation. Making this ae eq III-25 becomes; 
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where sensitivities are now positive numbers. 


If a controllable gain coefficient could be introduced into the numerator, 










so that an attenuation of the forced error could be provided, then the system 

might be made to perform without excessive errors. This could be done, it 
appears, if the wind acceleration term could be attenuated after entering the 

system. This cannot be done ; but an approximation of this desired result 

is achieved if additional feedback loops are provided as indicated in figure 

V-4. The performance equation for this system, when (PF) (aero) (ind) is 

assumed to be exactly equal to (PF) (aero)[AA ] 
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n the derivative auxiliary feedback sensitivity is set at zero, this becomes, 


h proportional plus first and second integral mechanization equation, 
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In this equation, a certain amount of gain control has been introduced 
o the numerator, at the expense of introducing a term giving a forced 
or caused by wind acceleration. Routh's stability criteria show that this 
stem, with proper selection of sensitivities, is stable. Limitations of 
ye did not permit the authors to make a complete analysis of this control 

. However, a sample result, selected at random is calculated. The 

ominator is first factored into two unequal quadratics, and the coefficients 
ared. They become: 
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These gave, using the equations given above, (DR), = 9.18, 
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he magnitude of the maximum forced error was found to be 28.5 nautical 
mile 3, and the disturbance which gave this maximum deviation was a sinu- 
bidally varying wind with a velocity magnitude of 265.5 nautical miles per 


our. This is a considerable improvement over the uncompensated system. 





Next, S(aux)[AA | was set at unity to remove the accleration term, and 
aux)[AA ] was set, at a random value. 


The performance equation, with quartic mechanization, becomes 
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Equating coefficients, as before, to those of a quartic with the desired 
characteristics, and selecting, at random, the following values: 


S(aux)[AA] = 70 


W, = Wo = 2W 


S oe 
Answers were obtained making (escicoutad = 0885 
CC(COVITA A] 


Seoncid 
(DOR), =(DR), =2.75, ccrcentsa | = 4.78 onal GIR = 3.25 


(C(CD)][AA] ae 
With the same wind input as that used in the preceding example, the maxi- 
mum error was found to be 48.8 miles. This, again, is an improvement 
over the uncompensated system. 

It appears, since in the first example the error was caused chiefly by 
the acceleration term, and in the second example by the second time rate 
of change of acceleration, that a combination of the two systems, varying 
both S(aux)(AA ] and S(aux)[AA | might give a satisfactory solution. Time 
prevents such a Study here. A simulator would greatly facilitate the work 
of selecting optimum sensitivities. 

B. Range Indication 

Lack of time prevented an elaborate examination of the response of 
the range indication system to large ranges of inputs. Since the track con- 
trol system indicates that the critical response is to sinusoidal inputs, a 
series of solutions was made for cubic and quartic responses to a sinu- 
soidal acceleration input with a period of thirty minutes. 

As indicated in the plots of Chapter IV, the response is improvedas the 
damping ratio of the quadratic in the cubic equation, or the two equal quad- 


_ratics in the quartic equation, is decreased. Ata damping ratio near 0.3, 
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4. Range Indication Mechanization Summary 


Examining the control equation obtained in Chapter IV, it is seen that 
it is possible to eliminate forcing function terms of orders up to and in- 
cluding wind acceleration terms. Ina stable system, higher order terms 
cannot be removed. A general performance function including all mecha- 
nization equations that eliminate acceleration and lower order terms can 


be written 


S(n=1)7 Ppt DS (n)T 


(PF) ech) (AA)(I AG) = i a 
re PSP + + Sino 


If more terms are included in the denominator, velocity and position 
errors result; if more terms are included in the numerator, control is lost 
over low order terms, while higher order terms are eliminated. If terms 
are omitted, instability results. 

A block diagram showing a method of mechanizing the performance 
function eq (x%~}9) follows, fig (V-5). 

For both track control and range indication systems,. the response for 
Sinusoidal forcing functions of 84.4 minute period appears to be the same 
regardless of the choice of damping or order of the equation in the 


mechanization system. 
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APPENDIX A 


DISCUSSION OF SIMPLIFYING ASSUMPTIONS 


Several simplifying assumptions were made in reaching the mechani- 
zation equations and performance equations discussed and analyzed in this 
thesis. These will now be examined. 

1. Omission of the Geodesic Acceleration Term 

The geodesic acceleration arises from the fact that navigation is 
performed on the surface of the earth, where the gravity equipotential 
forms the geoid. This geoid may be considered to be an ellipsoid with 
roughness. The deflection of the vertical due to this roughness never ex- 
ceeds seventy seconds of arc, and (especially over land masses) is usually 
far smaller. 

While the missile travels over this ellipsoid, computation assumes 
that the surface is a sphere. The shortest distance between two points on 
this sphere is a great circle. The gravity verticals fall in this great circle 
plane. The great circle plane transfers to the surface of the ellipsoid as 
a plane section, but the gravity verticals no longer lie in this plane, but in 
an arc intersecting it at the points of departure and destination. The mis- 
sile, traveling along this curved path, experiences a horizontal accelera- 
tion. The radius of curvature of this path, which is a minimum at 45 degrees 
north or south latitude, is never less than one million miles. The geodesic 
acceleration is therefore of negligible magnitude for missiles with veloci- 


ties of a few thousand miles an hour. 
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Two figures (Fig A-1 and 2) prepared by Dr. W. Wrigley of the 

M.1I.T. Instrumentation Laboratory, illustrate this geodesic acceleration. 
2. Compensation for Coriolis Acceleration 

The acceleration of the Coriolis arises from the fact that the mis- 
sile is moving in earth space, which possesses an angular velocity with 
respect to inertial space. This is demonstrated in the accompanying fig- 
ure (Fig A-3). Compensation is independent of the heading of the missile, 
depending only on a knowledge of the earth's angular velocity with respect 
to inertial space, which is accurately known, on latitude and on ground 
speed. An indication of latitude and of ground speed are available within | 
the missile, the accuracy of which are dependent on the success of the con- 
trol system. Errors resulting from inaccuracies in these measurements 
will be negligible in any practicable system. 

A certain amount of cross-coupling of the track control and range 
systems will occur through the Coriolis Computer, but it will remain small 
if the component of missile velocity perpendicular to the track is not large. 

3. Small Angle Pendulum Assumption 

For simplicity of calculation, a linearization assumption has been 

made in the output of the pendulum. The pendulous accelerometer actually 


solves the equation 


*[E -(cm)](hor) _ tan (DC) ¢rue) 
BIR 


The control equations are based on the assumption that 


a 
[E = (cm) J(hor) = (DC) tue) 
Em 
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The validity of this assumption has been tested, by comparing the out- 
put, for the track control system, to a pulse of wind acceleration of 2700 
feet/ min” , lasting for one minute, using the transcendental and the simpli- 
fied equation. The results are shown in the accompanying figure (Fig A-4). 

If necessary, a tangent pick-off can be used on the accelerometer in- 
stead of the linear pick-off postulated in this discussion. 

4. Neglect of Vertical Acceleration 

If the missile does not maintain absolutely constant altitude (or if 

the magnitude of the gravity force varies), the pendulum angle will not be 
a function of horizontal acceleration alone. This effect, together with the 
effect of an aerodynamic lag and a pendulum lag, is considered in the fol- 
lowing treatment of the range indication problem. There are changes in 
some of the geometric equations, as indicated: 


Aerodynamic lag term: 


Ore- (cm)) (hor) 1 (A-1) 
Vi@iry-E1 (CT aero) Ph 


Accelerometer lag term: 


(DC) truexcorr) (CT) pu) ee t 


Ope -(cm)I(hor) (A-3) 


RE 


r (true) 
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The effect of vertical acceleration becomes very clear. If the mis- 
sile momentarily loses lift and ''drops"’, the denominator goes to zero. 
The numerator, meanwhile, acquires a term which gives a forced error 
from wind acceleration. This error will change sign with change in the 
direction of vertical acceleration, but the effect of the denominator will 
keep the opposing effects from balancing. The damping in the denominator 
cubic, however, will return the missile to the track after the vertical ac- 
celerations have passed. 

Even with a downward acceleration greater than the effect of 
gravity, the system retains stability. 

Reasonable values for (CT) (26 ro) and (CT) pu) will vary from a 
fraction of a second to a few seconds. These will, therefore, have little 


effect on the system. 


5. Assumption of Perfect Aerodynamic Response 

Paragraph 4 indicates that neglect of aerodynamic response terms 
can have little effect on the performance of the longitudinal system. In 
the track control system, however, with a fourth order equation or higher, 
the system appears stable if a positive pendulum angle causes either a 
left or a right rudder deflection. It seems, therefore, that even a small 
aerodynamic lag might cause the system to become unstable. The problem 
was therefore recalculated for one choice of (DR), using, instead of perfect 
aerodynamic response, a pair of equations. The first gives the desired 
heading of the missile as a function of pendulum angle; the second relates 


this desired heading to the actual heading of the missile through an aero- 
dynamic lag. 
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The variation of the response of this system from the idealized 
system discussed in Chapters II and IV could not be detected in the plots 
drawn by the Differential Analyzer. 
6. Assumption of Perfect Gyros and Accelerometers 
Since this system is an angle measuring system which carries as 
a reference an indication of the point of departure dependent on the accu- ' 
racy of the gyros, any drift in the gyros which define the vertical plane 


will give an equal angular error in range and deflection. The third gyro, 


- A) If (DC)anarteyat dt 


(A-7) 
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* which maintains the orientation of the track plane, will introduce a non- 
linear error. The gyros have been assumed perfect primarily for sim- 
plification of the problem. If gyros of sufficient accuracy are found to 
be impracticable, the necessary inertial reference can be found from ob- 
servation of the fixed stars. 

Permanent errors in the zero indication of the pendulums also 
introduce an error on the earth's surface equal in — to the pendulum 
error. 

Single degree of freedom gyros and pendulous accelerometers of 
types which may prove satisfactory are being developed in the Instrumen- 
tation Laboratory at M.I.T. These units are shown diagrammatically in 
Figs (A-6) and (A-7). 

7. Omission of Constants of Integration 

Wherever possible, throughout this thesis, constants of integration 
have been assumed to be zero. For example, in the range indication prob- 
lem, the missile has been assumed to have no airspeed. This does not 
change the time response or the size of the range or track error resulting 
from wind disturbances, but will have a large effect on the distance 
travelled while the system is recovering from a large disturbance. 

At all times, the system has been assumed to be in equilibrium 
when hit with disturbing wind accelerations in the time solutions of Chap- 
ters III and IV. 
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APPENDIX C 


GLOSSARY 


The system of notation adopted for use in this thesis has been selected 
because it fulfills the following desirable objectives: 

1. It is easily learned 

2. It is adaptable to a wide range of situations 


3. It is built up almost exclusively of characters found on the 
keyboard of a standard American typewriter 


4. Any one of the compound symbols of the system is readily 
interpreted without recourse to an extensive glossary 


A few simple examples Suffice to explain the operation of the notation 
system. These examples are shown in table C-3. The short table of key 
symbols given at the end of this expinaea will then furnish sufficient 
information to enable any compund symbol to be correctly interpreted, and 
to provide. the necessary tools for the generation of new symbols. 

A representative list of the primary symbols is given in the following 
table. 
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a 

A 
(CT) 
[((CT)PR] 
(C) 
(DR) 
(Dir) 
(DC) 
(FR) 
g 
(PF) 
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Acceleration 

Angle 

Characteristic time 
Characteristic time - period ratio 
Correction 

Dampting ratio 
Direction 

Dynamic Correction 
Frequency ratio 
Gravity 

Performance function 
Radius 

Sensitivity 

Signal 

Time 

Velocity 

Vertical direction 


Frequency 


PRIMARY SYMBOLS 


TABLE C-l 
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(aero) 
(air) 
(app) 
(aux) 
(CD) 
(cm) 


(dep) 


(dest) | 


(hor) 
(ind) 
IR 

(kin) 


(long) 


NE 


CONF 


Aerodynamic 
Air mass 
Apparent value 
Auxiliary 
Control direction 
Controlled member 
Departure 
Destination 
Earth 

Horizontal 
Indicated value 
Inertial reaction 


Kinematical 


“T 
oo 


o 


Longitudinal plane; in 


range direction 
Missile 
Earth natural 


TIAL 


(pend) 
(pu) 


(ref) 
(res) 
(tc) 


(true) 


(Vert) 


Performance operator (d/dt) 
Pendulum 

Pendulous unit 

Range 

Reference 

Resultant 

Time 


Track Control; in track con- 
trol plane 


True value 
Vertical 


Axis through nose and c.g. 
of missile 


Axis through right wing and 
c.g. of missile 


Axis through c.g. perpendicu- 
lar to XY plane, directed 
downward 


MODIFYING SYMBOLS (SUBSCRIPTS) 


TABLE C-2 
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Symbol Definition 
Vig - (cm)] Velocity of the controlled member with 


respect to the earth 


S ; Sensitivity for correction-of the control 
[C(CD) ][AA] direction with angle input and rate of 
change of angle output 


(C) A Correction to the angle of the vertical 
(Vert)(tc) in the track control plane 
Wyre Natural earth frequency equal to 
g 
IR 
A Angle about the Z axis of the controlled 
Z(cm) member 
(DC) Indicated dynamic correction of track 
(ind) (tc) control system 
(DC) Indicated dynamic correction of longi- 
(ind) (long) tudinal system (i.e. range indication 
system) 
As (ind) Indicated range angle 


Typical Examples Showing How Symbols are Compounded from the Elemen- 
tary Forms 


TABLE C-3 
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